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ABSTRACT 

We consider the phase diagram of a classical fluid in the presence of a random pinning potential 
of arbitrary strength. Introducing replicas for averaging over the quenched disorder, we use 
the hypernetted chain approximation to calculate the correlations in the replicated liquid. The 
freezing transition of the liquid into a nearly crystalline state is studied using a density func- 
tional approach, and the liquid-to-glass transition is studied using a phenomenological replica 
symmetry breaking approach introduced by Mezard and Parisi. The hrst-order liquid-to-crystal 
transition is found to change to a continuous liquid-to-glass transition as the strength of the 
disorder is increased above a threshold value. 



PACS Numbers: 64.70.Pf, 64.70.Dv, 05.20.Jj 



The equilibrium phase diagram of a classical system of interacting particles in the presence 
of quenched pinning potential is a subject of much current interest, in view of understanding the 
behavior of various systems such as magnetic bubble arrays ffl, Wigner crystals of electrons Q, 
and flux lines in the mixed phase of high-Tc superconductors |y] . In particular, for layered type-II 
superconductors such as Bi2Sr2CaCu20s in a magnetic field H perpendicular to the layers, the 
(T, H) phase diagram is especially interesting. The flux lines in these materials may be viewed as 
columns of interacting "pancake" vortices residing on the layers, and the properties of the mixed 
phase may be described in terms of the classical statistical mechanics of these point-like objects. 
At low enough fields, a first-order melting transition separates an ordered vortex lattice from a 
disordered "vortex liquid" state Q . Existing theoretical studies |^ suggest that weak point 
disorder only slightly distorts the crystalline state, leading to the observed "Bragg glass" phase 
with quasi-long-range translational order The freezing line of a layered system of pancake 
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vortices in the presence of weak point pinning has been studied theoretically assuming the 
solid phase to be a harmonic crystal. When the magnetic field is increased, or equivalently with 
increasing disorder Q, the experimental transition becomes continuous |Q, |l^, and the "Bragg 
glass" is replaced by an amorphous state, which is believed to be a vortex glass |11, 12 1. 



This scenario could be a very general one, e.g. the first-order liquid-to-crystal transition 
in a three-dimensional (3d) fluid may be driven by quenched disorder into a continuous liquid- 
to-glass transition. In the present Letter, we consider, rather than vortices, a simpler and 
better documented system, namely a fluid of hard spheres in a random pinning potential of 
arbitrary strength. This is also a richer system since it offers the possibility of an intrinsic, 
though metastable, glassy phase, and the stabilization of this phase by quenched disorder is 
of considerable interest. We expect that our results could be generalized and applied to the 
systems listed above. Using two "mean- field"- type approaches based on the "replicated liquid 
formalism" ^ |l^ , we obtain a phase diagram in the density ( p) - disorder ( 6) plane which 
shows crystalline, liquid and glassy phases. One expects that : i) the transition between crystal 
and liquid phases remains first-order, while the transition point itself is shifted with increasing 
disorder; ii) from earlier work |^ and arguments |15] based on the Lindemann criterion, the liquid 
phase is favored by the disorder. Our phase diagram is consistent with these expectations. We 
also find that the first-order crystallization transition is replaced by a continuous glass transition 
as the disorder strength is increased above a threshold value. 

We consider a fluid of monodisperse hard spheres in the presence of an external, short-range, 
random pinning potential $(r). In the absence of disorder, the crystal becomes the stable phase 
at r/ = r^j ~ 0.49 {p ~ 0.94) where the packing fraction rj and the dimensionless density p 
are related by rj = 7rp/6 (spheres of diameter 1). If the hard-sphere fluid is kept in a "super- 
compressed" state, then slow dynamics effects, analogous to those occurring near the structural 
glass transition in supercooled liquids, are observed |16]. Whether a true thermodynamic glass 
transition occurs in this system at a packing fraction lower than that for random close packing 
{rjrcp — 0.64) is unclear |17]. In analogy with experiments, one may define the "glass transition" 
to occur at the point where the characteristic relaxation time in the fluid exceeds a given value 
Tc- Then, the phase diagram of the pure hard-sphere fluid would show a crystalline phase that 
is stable for ry > 0.49, and at r/ = r]g, just below rjrcp, a transition from the (metastable) liquid 
to a (metastable) glassy phase. 

The random potential <I>(r) is assumed to be gaussian with an exponentially decreasing 
correlator V{f — r') = $(r)<^(r') = A • exp[— (|r — r'p)/^^]; <I> = 0, (~ means averaging over 
the probability distribution of <I>), and a correlation length ^ = 0.25. This choice implies that 
local minima of ^ contain at most one sphere. The parameter A measures the strength of the 
disorder. As soon as an external potential is introduced, the temperature becomes a relevant 
parameter : the disorder may be considered perturbative ii 6 = /3^A <C 1; /3 = l/ksT, while 
(5 S> 1 corresponds to the strong pinning limit. 

We make use of the "replicated liquid formalism", introduced for studying fluids in porous 
media [^], flux-lattice melting and more recently, the structural glass transition |14]. The 
replica method is used for averaging over the disorder and mapping the problem of an inhomoge- 
neous liquid in a random potential onto the problem of a homogeneous, multicomponent mixture 
on which standard tools of liquid theory can be applied. Before averaging over the disorder, 
the particles are coupled to <I>(r) and interact via a hard-sphere repulsion. After introducing 
n copies ("replicas", labelled a, 6 = l,...,n) of the system and averaging over the gaussian 
disorder, one is left with a mixture of n components, interacting with v.t{r) = —(5V{r) plus a 
hard-sphere repulsion if the two particles belong to the same replica, and with VQ{r) = — /3V(r) 
otherwise. 

The structure of the replicated liquid is expressed in terms of the pair correlation functions 
9ab{f) = l + ^afel'') between particles of replicas a and h. Direct correlations functions Cab(r) are 
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introduced, related to the hat via the Ornstein-Zernike relation, expressed in Fourier space as 

Va, b; hab{q) = Cab{q) + ^ Kdiq) ■ Pd ■ Cdb{q)- (1) 

d 

Following Refs |, ||], we use the closure scheme knows as the hypernetted chain (HNC) ap- 
proximation in which the n — > limit can be taken analytically. In our case, the HNC equations 
read : 

gabir) = (1 - 6ab + dab ' e(r - 1)) • exp[/32v(r) + hab{r) - Cab{r)]. (2) 

At relatively low values of the packing fraction, replica symmetry is a natural assumption. Then, 
the set of functions gab {hab and Cab) reduces to g'* {h^, and c=k) \i a = h and go (^o and cq) if 
a ^ b, while all the densities pa are equal to p. The function g*(r) looks very similar to the 
usual pair correlation function of the pure system. In the absence of disorder, (7o(^) = 1) and it 
starts showing structure when b is increased. It exhibits peaks at the maxima of gi:{r), plus an 
extra peak at the origin induced by the attractive coupling between replicas. 

The freezing into the nearly crystalline Bragg glass phase occurs in the replica symmetric 
(RS) regime. According to the standard approach of density functional theory |^, the 
crystal manifests itself through a (almost) periodic modulation of the density Pcri^, obeying 
the approximate self-consistency equation: 



Pcrir) = p-exp 



dfCeif- r') ■ {Pcr{r) - p) 



(3) 



where the function Ce{r) stands for Y^b '^ab{^) = c* + (n — 1) • Co(r) = c^{f) — Co(r) for n — > in 
the RS regime. The determination of the freezing line in the presence of disorder then reduces 
to a standard density functional calculation with a modified direct correlation function. 

The glass transition in this system is studied using the phenomenological approach of Mezard 
and Parisi (MP) who found the occurrence of one-step replica symmetry breaking (RSB) 
in equations (|l|),(^ in the absence of external disorder and interpreted it as the emergence of a 
glassy phase. MP consider a functional of the pair correlation functions gab{f) which yields 
the set of equations (|l]),(0) upon functional differentiation with respect to Qab- The one-step 



RSB solution consists, as usual, in forming n/m groups, each containing m replicas [19|. Pair 
correlations Qab are set to be g^, for a = 6, gi for a 7^ 6 with a, h in the same group, and go if a 
and b are in different groups. Then, one looks for solutions making the free energy !F stationary 
with respect to g*,gi,go and m, with < m < 1 and n = 0. By analogy with one-step RSB 
in mean-field spin-glass models with multispin interactions, a "dynamical transition" density 
Pdyn is defined as the minimal density for which a RSB solution with m = 1, stationary with 
respect to g*,gi,go but not with respect to m, exists. The occurrence of a stationary solution 
with m = 1 at a higher value of p signals a thermodynamic glass transition. 

The approach of MP relies strongly on the assumption that in the vicinity of the glass 
transition, the configuration space is split into an exponentially large number of "metastable 
states" [^]. These metastable states come into existence at the "dynamical transition" density 
Pdyn, and the thermodynamic glass transition is expected to occur at the density where the 
configurational entropy associated with the metastable states vanishes pO| , ^] . Since the concept 
of metastable states is itself a mean-field one borrowed from the study of infinite-range spin- 
glass models, it is not clear whether this description survives in real 3d fluids. Nevertheless, we 
believe that both thermodynamics and dynamics of realistic 3d fluids near the glass transition are 
dominated by long-lived configurations [23, ^|, which would be the 3d analog of the metastable 
states. Then, the method proposed by MP may be regarded as a phenomenological way of 
finding when these long-lived states start playing a significant role. We also note that results 



very similar to those of MP have been obtained in a calculation |22] that uses a different method 
for locating the glass transition in the pure hard-sphere system. 
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Due to the finite dimensional character of the system, the "overlaps" gi{r), go{r) have a 
spatial structure which accounts for the short-range ordering of the particles in the metastable 
states. Roughly speaking, gi{r) describes short-range correlations of the average local density in 
a typical metastable state, and go{r) represents such correlations between different metastable 
states [go{r) = 1 in the absence of disorder]. We define the scalar Q = [/ dr 47rr^ {giii") " 
g(o(r))^]^/^ as a global order parameter. At zero disorder, the RSB transition occurs at pg ~ 1.16 
and is discontinuous, with a jump in the value of Q at the transition |]25| . 



We have solved the replicated liquid-state equations using the method proposed by Zerah [ 26 1 , 
and a grid of 512 points. Fig. 1(a) presents our phase diagram. The thermodynamic glass 
transition found in |14| extends (thick line) from p ~ 1.16, 5 = to a "tricritical" point T : 
pt ^ 1.11, 6t ^ 1.2. For 6 < 5t, Q jumps discontinuously at the transition, while for 6 > 5t (thin 
line), Q grows continuously from zero as p is increased across the transition line. The dashed line 
represents the freezing line obtained from the density functional calculation. The liquid phase 
is favored by the disorder and the freezing line crosses the glass transition line at C: pc — 1.10, 
6c ~ 1.67. The glass phase is also favored by the disorder until 6 ~ 1.7 where the liquid shows a 
re-entrance. The upper part of the freezing line lies in a region where the disorder is not weak 
{5 > 1) and so, the stability of the Bragg glass phase is not ensured. However, if we assume that 
the ordered phase remains stable in this region, then extrapolation of the freezing line into the 
glassy domain leads to the inset of Fig. 1(a), where only stable phases are shown. As the density 
is increased at constant 6, the system undergoes a first-order transition to a nearly crystalline 
state if (5 < 6c, and a continuous glass transition for 6 > 6c- Thus, the phase diagram exhibits a 
multicritical point where a line of continuous liquid-glass transition ends at a line of first-order 
transitions representing a liquid-crystal transition for 6 < 6c and a glass-crystal transition for 
6 > 6c. 

We now provide some of the technical details of our calculations. The HNC free-energy per 
unit volume V reads : 
2/3.F f 

—— = {g^ (In g^ - 1 + (3v^) + {m - I) gi {In gi - 1 + Pvo) - m go {In go - 1 + (3vo)} 

nV J 

I pK - ^{hl + (m - l)hl - mhl) + — — ln(l + pK - phi) 

-— ln(l + pK + (m - 1) phi - m pho) - — — - , , ~- \ ■ (4) 

In the n ^ limit, a physically stable solution must be a minimum with respect to variations of 
g^f , but a maximum with respect to gi,go- As m — > 1 , the equations for gi decouple from those of 
g^, and go, and one can get simultaneously the RSB solution {g>t,gi ^ go) and the RS one {g*,gi = 
go). Both stable RS and RSB solutions exist in the domain bounded by the lines (IN,DT, 6 = 0) 
in Fig. 1(b). The thermodynamic glass transition (line TGT) occurs when dJ- / dm\m=i = 0, i.e. 
when the second term (^'[g^,,gi,go\ of the expansion = (t[g^,, go\ + {m — 1) ■ + 0{{m — 1)"^) . . . 
vanishes (this is our practical criterion for determining the location of the line TGT). 

Examples of RS and RSB solutions on two sides of the line TGT aX, 6 = 0.3 are shown in 
Fig. 2(a). These plots illustrate the discontinuity in Q across the line TGT. This discontinuity 
decreases with increasing 6, and the transition eventually becomes continuous for 6 > 6t on the 
line CT (continuous transition). Fig. 2(b) presents the RSB solution next to the line CT at 
6 = 1.8. The functions gi and go look very similar, indicating that Q vanishes at the transition, 
as shown in the inset. The transition ai 6 > 6t occurs via a bifurcation mechanism. This can be 
understood by expanding to second order in \g^,A.gi, A^o around the RS solution. Defining 
Afci = \/m{l — m)(A^i — A/iq), A/cq = (1 — m)A/ii + mA/iQ, and A = 1 + ph^ — pho, we get : 

/3A^ 



nV 



J{Ah„Ako) - />V4 11 df {g^' - l)Akl + J dq/{2n)^ A'^Ahf^ , (5) 
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where ^7 is a m-independent quadratic form. We have diagonahzed a discrete version of the 
quadratic form enclosed in the braces, and found that its lowest eigenvalue Xmin changes its 
sign, from positive (in the liquid phase) to negative (in the glass phase), just when Q vanishes 
(inset of Fig. 2(b)). The RS solution turns unstable as Xmin becomes negative, providing a 
precise determination of the line CT. This line admits a continuation for 5 < 6t, labelled IN 
(instability line) in Fig. 1(b), which does not seem to have any physical meaning. A similar 
computation around the RSB solution enables us to determine precisely the bifurcation point 
at which the RSB solution, m = 1, appears. This calculation yields the "dynamical transition" 
line DT. As the transition becomes continuous near 5 ^ 5t, finding numerically the line TGT 
becomes very difficult because the ^' variations become vanishingly small. However, the fact 
that the lines CT, DT and IN converge to the same point T leads us to the conclusion that the 
TGT line also ends at T, as depicted in Fig. 1(a). 

The pure hard-sphere fluid freezes into a fee lattice. We use the Ramakrishnan-Yussouff 
method [|l^] for expressing the density pcr ('^) in terms of three "order parameters" ~ the Fourier 
components, /U3,i,i, of the density for two sets ((1, 1, 1) and (3, 1, 1)) of reciprocal lattice 

vectors, and r], the fractional change in the average density: 



In^ 
P 



VP Ce{0) + Pl,l,l ^e(|i^l,l,l|)e''^'''''" + /^3,l,l E Ce{\K3,l,l\y'''-'-'- 

(1,1,1) (3,1,1) 



(6) 



We believe that this parametrization is justified in the presence of weak disorder because the 
nearly crystalline Bragg glass phase exhibits well-defined peaks in its structure factor. The HNC 
direct correlation function is used as input, and the coexistence line is obtained from the usual 
thermodynamic criterion ||T^: 



J dr i^p^^.ln{P^)-p,r+p^-lJdrJ dr' Ce(f- f ') • (P^r " P) ' (Pc. - p) = 0. (7) 

At (5 = 0, we get pf = 0.93, with rj ~ 15%. As 5 is increased, pj increases (see line FL in 
Fig. 1(b)) and rj decreases, reaching a value of about 6% near the point C. 

To summarize, we have used a combination of the replica method, liquid theory and density 
functional theory to obtain the phase diagram of a simple classical fluid in a random pinning 
potential. Our calculations quantify the effects of the disorder on the crystallization transition of 
the pure fluid and its glass transition in the metastable "supercompressed" regime. We find that 
the first-order crystallization transition of the pure fluid changes to a continuous glass transition 
as the strength of the disorder is increased above a critical value. The phase diagram we 
have obtained looks qualitatively similar to that of layered type-II superconductors if, instead of 
increasing the density p, one decreases the temperature, and if one replaces the disorder strength 
6 by the magnetic field H. Since our calculations are mean-field in nature, they do not provide a 
conclusive answer to the question of whether a thermodynamic glass phase exists in 3d randomly 
pinned classical systems. Further investigations of this issue and extensions of our calculation 
to layered superconductors and other physical systems would be very interesting. 
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Figure captions 



Fig. 1: (a) The phase diagram in the density p - disorder 5 plane. T is a tricritical point where 
the nature of the glass transition changes from first-order to continuous, and C is a critical 
endpoint where the continuous glass transition line meets the first-order freezing line. The inset 
depicts the phase diagram in which only the thermodynamically stable phases are shown, 
(b) Various transition lines (see text for details) in the p — 6 plane. TGT: first-order glass 
transition for weak disorder; CT: continuous glass transition for strong disorder; DT: dynamical 
transition; IN: instability of the RS solution; FL: freezing line obtained from density functional 
theory. 

Fig. 2: (a): HNC pair correlation functions near TGT, at 5 = 0.3. Thick and dashed lines stand 
respectively for gi and go at p = 1.15, slightly on the glass side, while the "diamond curve" is the 
replica symmetric solution at p = 1.14, on the liquid side. Note the jump of Q at the transition. 
Inset : g^, gi, go at 5 = 0.3, p = 1.15. 

(b): Pair correlation functions near the continuous transition CT, at 5 = 1.8. Inset : The "order 
parameter" Q (see text) which vanishes at the transition point. Also shown is A, the lowest 
eigenvalue of the quadratic form enclosed in braces in Eq.(^, which changes its sign at the 
same point. 
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(b) 

Fig. 1 
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Fig. 2 



